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Recent status of neutrino oscillation phenomenology with four neutrinos is re- 
viewed. It is emphasized that the so-called (2-|-2)-scheme as well as the (3+1)- 
schemc arc still consistent with the recent solar and atmospheric neutrino data. 



1 Introduction 



There have been several experinientsldl^'BHi3^BEIi3^Eti]'Eilli3£j which suggest neu- 
trino osciUations. To explain the solar, atmospheric and LSND data within 
the framework of neutrino oscillations, it is necessary to have at least four 
kinds of neutrinos. It has been shown in the two flavor framework that the 
solar neutrino deficit can be explained by neutrino oscillation with the sets 
of parameters {Am%,sin^ 29q) ~ (^(lO^^eV^), ©(lO^^)) (g^^ (smaU mix- 
ing angle) MSW solution), (O(lO-^eV^), 0(1)) (LMA (large mixing angle) 
MSW solution), (C'(10-W),O(l)) (LOW solution) or (©(IQ-iOeV^), 0(1)) 
(VO (vacuum oscillation) solution). At the Neutrino 2000 Conference the Su- 
perkamiokande group has updated their data of the solar neutrinos ^and they 
reported that the LMA MSW solution gives the best fit to the dataEI. At the 
same time they also showed that the scenario of pure sterile neutrino oscilla- 
tions i>e is excluded at 95%CL. It has been known that the atmospheric 
neutrino anomaly can be accounted for by dominant <-> Vt oscillations with 
almost maximal mixing (Am^^.j^, sin^ 2fiatm ) ~ (10-2-5eV^L0). Again the 
Superkamiokande group has announced 113 that the possibility of pure sterile 
neutrino oscillations <-> i>s is excluded at 99%CL. On the otheLhand, com- 
bining the_final result of LSND ancLthe negative results by E776tj {vp, Vg), 
Karmen2E3 (t'p — > Ve) and BugeyEJ {ve v^), the oscillation parameter sat- 
isfies 0.1 eV^ < Arn^ <8 eV^ and 8 x 10-'*< sin2 26'<0.04 at 99%CL. In this 
talk I will review the present status of four neutrino scenarios in the light of 
the recent Superkamiokande data of the solar and atmospheric neutrinos. 



2 Mass patterns 

In the case of four neutrino schemes there are two distinct types of mass pat- 
terns. One is the so-called (2-|-2)-scheme (Fig. 0(a)) and the other is the 
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Figure 1: Mass patterns of four neutrino schemes, (a) corresponds to (2+2)-scheme, where 
either (lAm^^l = Am|, lAmj^l = Am^^^^) or {\Aml^\ = Am|,, \Aml^\ = Am^^^). 



(b) and (c) arc (3+l)-scheme, where |Am,|-|^| = Am^gj^^ and either (lArrtj-^l = Am^ 



I Am 



ArriQ, 



-411 
I Am? 



^'"atm 



0' 



) is satisfied. 



(3+l)-scheme (Fig. |^(b) or (c)). Depending on the type of the two schemes, 
phenomenology is different. 



2.1 (3+l)-scheme 

It has been shown in Refs. 0E1 using older data oIXSNdEI that the (3+1)- 
scheme is inconsisteiLL with the Bugey reactor datalj and the CDHSW disap- 
pearance experiment^ of v^. Let me briefly give this argument in Refs.Hilj. 
Without loss of generality I assume that one distinct mass eigenstate is (See 
Fig. 11(b) or (c)) and the largest mass squared difference is ATO43 = ^"^lsnd- 
In the case of (3+l)-scheme the constraints from the Bugey and CDHSW 
data are given by 

1 - P(Pe ^ Ve) = 4|C/e4p(l - \Uei\'')^A:i < siu^ 20Bugcy ( Am^g ) A43 , 

1 - P{iy^ ^ ly^) = 4|C/^4|'(1 - \Uf,4\^)A43 < sin^ 20cDHSw(Am23)A43, 



respectively, where A43 = sin^ (Am|3L/4i?), sin'^ 26'Bugoy and sin^ 20cdhsw 
stand for the values of the boundary of the excluded region in the two flavor 
analysis as functions of Am^ (See Fig. |^). To explain the solar neutrino 
deficit and the zenith angle dependence of the atmospheric neutrino data it is 
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Figure 2: The allowed region of LSND from the final result (the inside of the thick and thin 
solid lines corresponds to the allowed region at 99%CL and 90%CL, respectively) and the 
excluded regions of E776, Karmen2, Bugey, CDHSW (the right hand side of each line). The 
right hand side of the line (Bugey+CDHSW) stands for the excluded region in the case of 

the (3+l)-schcmc. Only the four isolated shadowed areas at Am^gj^^ ~0.3, 0.9, 1.7, 6.0 
eV'^ are consistent with the LSND allowed region at 99%CL in the (3+l)-scheme. 
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necessary to have |?7e4p < 1/2 and |C/^4p < 1/2 and therefore I get 



\UeA\' < 



1- Jl- sin^ 2e'Bugcy(Am23) 



1 - Jl - sm^ 26'cDHSw(Am|3) 



(1) 



On the other hand, the appearance probabihty P{Vf_i — ^ I'e) of LSND in our 
scenario is given by 



P{i^^ ^ Ve) = 4|t/e4nC/MpA43 = sin^ 26'lsnd(Ato|3)A43, 



(2) 



where sin^ 20LSND(Am43) stands for the value of sin^ 29 in the LSND aUowed 
region in the two flavor framework. From (|l|) and @) I obtain 



sin^ 26'LSND(Am43) < 



I - sin^ 26'Bugcy(Am|3) 



1 - a/ 1 - sin 26'cDHSw(Am|3) 



(3) 



The value of the right haiid side of (||) is pfotted in Fig. g together with the 
aUowed region of LSND Ej. At 90%CL the aUowed region of LSND does not 
satisfy the condition ^ for the_(3+l)-scheme, an|d,actuaUy it used to be the 
case with older data of LSND Il3 even at 99%CLll§lia. However, in the final 
result the allowed region has shifted to the lower value of sin^ 29 and it was 
shownEy that there are four isolated regions AmLg^D —0-3, 0.9, 1.7, 6.0 eV^ 
which satisfy the condition 



2.2 (2+2)-scheme 

In the case of the (2+2)-scheme, assuming the mass pattern in Fig. |l| (a) with 
Amli = Atoq, Am^2 = Am^^j^, ATO43 = ^^lsnb^ the constraints from the 
LSND, Bugey and CDHSW data are given by 

Pii^^ ~* i>e) = MUesU;^ + UeiU;^\^A32 = siu^ 20LSND ( Am22 ) A32 , 
1 - Pii?e ^ J>e) = 4(|C/e3p + \Uei\''){l - |C/e3p - \Uei\'')A32 

< Sin^ 26lBugcy(Am22)A32, 

1 - P(i.^ ^ i.^) = 4(|C/^3p + \U^4\^){1 - |f/^3|' - iC/^nAsz 

< sin2 26icDHSw(Am^2)A32, (4) 

where A32 = sin^ (Ato^2-^/4£^)- 
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Figure 3: The allowed region of the (2+2)-scheme with i5i = 7r/2 at 90 %CL constrained 
by the Superkamiokande atmospheric neutrino data. The shadowed regions stand for the 
allowed regions from the combined analysis of the solar and atmospheric neutrino data. 
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It has been shown HEj that these conditions are consistent with aU the 
constraints from the accelerator, reactor data as well as solar and atmospheric 
neutrino observations. As I will show, to account for both the solar neutrino 
deficit and the atmospheric neutrino anomaly, it is necessary to have 

m,3\' + - |C/^3p - \U,,\') - Oil), (5) 

so I take Am§2 as small as possible, i.e., Am§2 = 0.3 eV^ so that be 
consistent with the CDHSW constraint (Q). 



3 (2+2)-scheme 

3.1 Analysis of the solar neutrino data 

The solp: neutrino data were analyzed in the framework of the (2+2)-scheme 
by Ref. E3 on the assumption that 11^3 = C/e4 — 0, which is justified from the 
Bugey constraint | ?7e3 P + 1 t^e4 P *C 1, and Am^^j^, AmLgN^ oo which is also 
justified since \Am1^j^/2E\, \Am^gj^j~,/2E\ :s> \/2GFNe for the solar neutrino 
problem, where Gp and stand for the_Eermi constant and the electron 
density in the Sun. The conclusion of Ref. E3is that the SMA MSW solution 
exists for < Cs <0.8, while the LMA MSW and LOW solutions survive only 
for < Cs <0.4 and < <0.2, respectively, where Cs = |f/siP + |C^s2p- 



3.2 Analysis of the atmospheric neutrino data 

The atmospheiic_neutrino data were analyzed by Jiefs. BE with the (2+2)- 
scheme. Refs.tMj assumed Ues = C/e4 = as in Ref.E3, and Atuq = was also 
assumed. Ref. E3 assumed Attilsnd — ^-^ so that the result with large 
|C/,i3p + \U^,4\'^ do not contradict with the CDHSW constraint (|). Ref. 
did not take into account the contribution fijajn AmLgND oscillation 
probability and their result is a subset-pf Ref.E3. 

Here I adopt the notation in Ref. for the 4x4 MNS matrix: 

UmNS = i?34(| -034)i?24(^24)i?23(|)C/23(023,^l)f/l4(^14,<53)f/l3(^13,'52)i?12(012) 

(6) 

where C/23 (023, 5i) = e^'^^^-^ R23{-923)e-^'^'^' , (714(^14,^3) = e^'^^^'''^^Ri4{9i4) 
g-V6.53Aia/2^ C/i3(0i3, <52) = e2»^^^«/^i?i3(0i3)e~2^'^^«/^, RjkiO) = exp (tTjkO) , 
(Tjk)em = * i^je^km - Sj„,5m), 2X3 = diag(l, -1, 0, 0), 2V3A8 = diag(l, 1, -2, 0), 
2V6A15 = diag(l, 1, 1, —3) are 4x4 matrices {Xj are elements of the su(4) gen- 
erators). With the assumptions Am^i = 0, |C/e4p = S14 = 0, |C/e3p = cf^sl^ = 
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0, Oi2, 013 and ^14 disappear from U and v^. decouples from other three neutri- 
nos. Thus the problem is reduced to the three flavor neutrino analysis among 
v^, Vr, i^s and the reduced MNS matrix is 

U= Ur2Ur3Uri ^ e'( f -«34) A7 £,-lg«e24 A5 g'C'^^ - f ) A2 ^ 
V Us2 Us3 Usi I 

with D = diag (e^''^/^, 1, e-''^^/^) (Aj are the 3 x 3 Gell-Mann matrices) is the 
reduced 3x3 MNS matrix. This MNS matrix U is obtained by substitu- 
tion 6*12 6*23 — 7r/2, dirh^ 0-24, di2 7r/2 — 6*34, 5 ^ 5i in the standard 
parametrization in Ref. 123. It turns out that ^34 corresponds to the mix- 
ing of Vfj, ^ Vt and Vy^ <^ Vs, while ^23 is the mixing of the contribution 
of sin^ (Am^tmi/4ii^) and sin^ (Atolsnd^/^-E') in the oscillation probability. 
The allowed region at 90%CL of the atmospheric neutrino data is given by the 
area bounded by thin solid linesJji Fig. ^ for 61 = tt/2. The allowed regions 
for Si — 0, 7r/4 are given in Ref. c^. 

3.3 Combined analysis of the solar and atmospheric neutrino data 

In Fig. I, the lines given by = \Usi\'^ + \Us2\'^ = IC23C34 + S23S34S24e*''ip = 
0.2, 0.4, 0.8 are depicted together with the allowjed region o£-|the atmospheric 
neutrino data. By combining the analyses of Ref.E3 and Ref.c3, I obtain the re- 
gion which satisfies the constraints of the solar and atmospheric neutrino data. 
The darkest, medium and lightest shadowed areas stand for I'atm + J^0(SMA, 
LMA or LOW), j/atm + i^0(SMA or LMA), t/atm + i^0(SMA), respectively. 
Although this result is not quantitative, it gives us a sense on how likely the 
(2-|-2)-scheme is allowed by combining the solar and atmospheric neutrino data. 
Let me emphasize that non-zero contribution of sin^ (AmLgNi3i/4_B) (i.e., the 
case of 6*23 > 0) to the oscillation probability is important particularly for the 
LMA ajid LOW solar solutions. The region of 623 > has not been analyzed 
by Ref.EJ. Let me also stress that both the solar neutrinos and the atmospheric 
neutrinos are accounted for by hybrid of active and sterile oscillations in the 
(2-|-2)-scheme. 

4 (3+l)-scheme 

After the work of Barger et al. §, people ii have investigated various conse- 
quences of the (3+l)-scheme. Here let me make two comments on the (3+1)- 
scheme. 
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4-1 Atmospheric neutrinos 

As in the case of the (2+2)-schenie, I assume Ue3 = Ud — Atoq — for 
simphcity. Then once again decouples from v^.^ v^, Vt and the probabiUty 
in vacuum for the atmospheric neutrino scale is given by 

^ ^t^) = Mu^3\Hi - \u^3\^ - \u^4\^)A32 + 2|t/p4p(i - Ic/mH: 

P{l^^ ^ l^r) = 43? [U^3U:^iU;^Ur3 + U;^UrA)] A32 + 2 1 (7^4 H C/r4 | ' , 

^ ^s) = 43? [U^3U:s{U;:,Us3 + U;^U,i)] A32 + 2|t/^4p|C/.4|', 

(7) 

where I have taken Am|2 = Am^^j^, Am43 = ^™lsnd ^.nd I have averaged 
over rapid oscillations: sin^ (AmLSND-^/4-£') ^ 1/2- Since the (3+l)-scheme is 
allowed only for four discrete values of Attilsnd) ™^ discuss Amf^g^^=0.3 
(|«7mP>0.34) and Amlsj,^=0.9 eV^ (|C/^4p ^ 0.03), 1.7 eV^ (|C/^4p ^ 
0.01), 6.0 eV^ (|t^Ai4p — 0.02), separately. For simplicity I assume Si=0 since 
the existence of the CP phase does not change the situation very much. 



4.1.1 Am2gj^j3=0.3 eV2 

Since we know from the Superkamiokande atmospheric neutrino data that the 
coefficient of sin^ (Am|^j^L/4i?) in v^) has to be large to have a good 

fit, I optimize 4|;7^3|2(1 - |J7^3|2 _ |c/^4|2) with respect to 6*23 for \U^,i\'^ = 0.34. 
When C/e4 = I have |[/^3|' ^023024, |t7^4p = 5^4, \UrA? = cl^cl^, |[/,4p = 
c\^s\^ in the notation of Ref.ll3, and it is easy to see 

4|;7^3|'(1 - |(7^3|' - It/^rfP) = 4sin2 2023 < 4 = 0.44, 

where equality holds when 6*23 — 7r/4. This is the value of ^23 for which the 
fit of the (3+l)-scheme to the atmospheric neutrino data is supposed to be 
the best for AmLgND^O.3 eV^. When 6*23 = 7r/4 the probability in vacuum 
becomes 

P{v^ -> Vr) = ^C24S34 - ^ Sin^ 2024C34^ A32 + ^£34 siu^ 2^24 

P{v^, Vs) = {cl^cl^ - ^ sin^ 26*24534^ A32 + ^534 sin^ 26'24. 

(8) 



From (H) 034 turns out to be the mixing of Vt and Vs as in the 

(2+2)-scheme. 
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I found from the explicit numerical calculationL3 that the fit of the (3+1)- 
scheme with eV^, |L'/j4p — 0.34, 6*23 — 7r/4 to the atmospheric 

neutrino data is very bad for any value of 634 and the region of Aml^g^^=0 .3 
eV^ is excluded at 6.9crCL. 



4.1.2 Am^sND =0.9, 1.7, 6.0 eV^ 

In this case |J7^4p<0.03 and I can put /7^4 = as a good approximation. 
Then the constant part in the oscillation probability disappears and this case 
is reduced to the analysis in the (2+2)-scheme with 623 = 0. The allowed 
region at 90%CL is given roug hly by -7r/4 < 6*34 < 7r/4, 0.8 < sin^ 26*24 < 1, 
where ^34 and ^24 stand for the mixing of <-> v^- and <-> Vs and the 
mixing of atmospheric neutrino oscillations, respectively. 



4-2 Oscillations of high energy neutrinos in matter 



When |J7e4p, |t^M4p ^^'^ |C^T4p are all small, it is naively difficult to distinguish 
the (3+l)-scheme from the ordinary three flavor scenario. However, because of 
the existence of the small mixing angles in Ue4, U^i and the large mass squared 



difference Ar 



'LSND 



the oscillation probability in matter can have enhancement 



which never happens in the three fiavor case. By taking 612 = 7r/4, 613 
023 = tt/A, 014 = e i\e\ « 1), 024 - S {\S\ « 1), 634 - 7r/2 in {§ I get 
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which is the same as the MNS matrix in Ref . up to the phase of each factor 



0, 



The probability P{i^f_, 



turns out to receive significant deviation from the 



vacuum one due to the matter effect for ~ C'(l) TeV, and the behaviors 



of 1 — P{i'fj. Vfi) are shown in Fig. IJ, where three cases of |ATOLgND|=0.9 
eV2 (8.8° < e < 12.2°, 6.4° <5< 8.9°), 1.7 eV^ (7.5° < e < 10.2°, 5.6° < 
5 < 7.7°), 6.0 eV2 (7.5° < e < 7.7°, 10.0° <S< 10.2°) are considered.^ The 
appearance channel which is enhanced is dominantly — )■ Vs, so it may be 
difficult to detect signs of this enhancement from observations of high energy 



" The eigenvalues of the mass matrix in this case turn out to be roots of a cubic equation and 
piialytic treatmeiitnpfithe oscillation probability is difficult, unlike the cases of three flavors 
CJ or four flavors EjtJi where one mass scale is dominant and the eigenvalues are roots of a 
quadratic equation. 
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cosmic neutrinos of energy E^, ~ 0(1) TeV, although this enhancement may 
be observed through neutral current interactions in the futurefl 




Figure 4: Enhancement of neutrino oscillations due to matter effect in the case of the (3+1)- 
scheme, where 823 = 45° is taken for all the cases. For each value of AmP= -0.9 eV^ 
(e = 10.0°, S = 7.5°), -1.7 eV^ (e = 5.6°, <5 = 7.5°), -6.0 eV^ (e = 6.4°, 5 = 8.8°), ten curves 
correspond to cos = —1.0, —0.9, ■ • ■ , —0.1 from the left to the right, where the zenith angle 
is related to the baseline L by L = — 2Hcos0 with i?=6378km. Most of the channel is 



5 Big Bang Nucleosynthesis 

It has been shownEl in the two flavor framework that if sterile neutrino have 
oscillations with active ones and if Am^ sin^ 26* > 3 x 10^'^ eV^ is satisfied 
then sterile neutrinos would have been in thermal equilibrium and the number 
N,y of light neutrinos in Big Bang Nucleosynthesis £JBBN) would be 4. This 
argument was generalized to the four neutrino caseEJ'LlI and by imposing all 
the constraints from accelerators, reactors, solar neutrinos and atmospheric 
neutrinos as well as the BBN constraint < 4.0 it was concluded that the only 
consistent four neutrino scenario is the (2+2)-scheme with the MNS mixing 

''Similar enhancement has been discussed-in a different context by Ref. 0. I thank Athar 
Husain for bringing my attention to Ref.Eil. 
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matrix 



Umns 



I Uel Ue2 Ue3 Uei \ / Cq Sq 

Uf_,i U^2 U^3 U^4: ^ e e 

Url Ur2 Ut3 Uri £ £ 

\Usl Us2 Us3 Usi J \~Sq Cq 



1/V2 I/V2 



(9) 



where Cq = cos 00, Sq = sinfl© and 9q stands for the mixing angle of the 
SMA MSW solar solution. In this case Cs = |C^siP + |C^s2p — 1 and the solar 
neutrino deficit would be accounted for by sterile neutrino oscillations ^ J^s 
with the SMA MSW solution while the atmospheric neutrino anomaly would 
be by <-> oscillations. This scenario is obviously inconsistent with the 
recent solar neutrino data by the Superkamiokande group, and the argument 
which has lead to (j^) has to be given up. 

Fortunately the upper bound on has become less stringent nam and 
Ni, — 4.0 seems to be allowed. Furthermore, it has been shown recently tn that 
the combinecLjanalysis of BBN and the recent data by BOOMERanGo and 
MAXIMA- IE3 of the Cosmic Microwave Background prefers higher value of 
N,y: 4 < A^^ < 13. Therefore all the four neutrino schemes of type (2+2) and 
(3+1) seemed to be consistent with the BBN constraint. 

On the other hand, it has been pointed out0 in the two flavor framework 
that for a certain range of the oscillation parameters neutrino oscillations them- 
selves create asymmetry between v and D and this asymmetry prevents Vg from 
oscillating into active neutrinos. Although this analysis has not been general- 
ized to the four neutrino cases, even if the upper bound of N^, becomes less than 
4.0 in the future, it might be still possible to have four neutrino schemes which 
are consistent with the BBN constraint as well as the solar and atmospheric 
neutrino data due to possible asymmetry in v and D. 



6 Conclusions 

In this talk I have shown that there are still four neutrino scenarios ((2+2)- 
as well as (3+1)- schemes) which are consistent with all the experiments and 
the observations, despite the recent claims by the Superkamiokande group that 
pure sterile oscillations i>e i^s in the solar neutrinos and pure sterile oscilla- 
tions <-^- Vs in the atmospheric neutrinos are disfavored. In particular, the 
reason that the (2+2)-scheme is consistent with the recent Superkamiokande 
data is because both solar and atmospheric neutrinos have hybrid oscillations 
of active and sterile oscillations. 
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